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HIGHLIGHTS 


► Two-fluid model with kinetic theory of granular flow is adopted for simulation. 

► Effect of gas flow rates on the multiple-spouted bed hydrodynamics is studied. 

► Effect of different particle sizes is investigated. 

► Spout diameter and fountain height increases initially and then decreases. 

► Static bed thickness has an obvious influence on the flow behavior in the bed. 
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In present work, the gas-solid two-phase flow in the multiple-spouted bed is simulated by the 
Eulerian-Eulerian approach, the Gidaspow drag model is chosen to describe the interface momentum 
exchange. The simulated results agree well with the experimental results from the literature. Then the 
detailed flow behavior in the multiple-spouted bed is investigated for 1.0 mm, 1.4 mm and 1.8 mm 
glass beads. The effect of the ratio of central/auxiliary and auxiliary/central gas flow on the 
hydrodynamic characteristics in the multiple-spouted bed is studied. The high velocity spouting gas 
dominates the flow pattern in the bed and restricts the development of the low velocity spouting gas. 
The distribution of voidage, the profile of particle velocity and the variation of particle volume 
concentration are obtained. The fountain height and spout diameter increases with increasing gas flow 
rate initially and then decreases. The increasing bed thickness has an obvious influence on the 
hydrodynamics in the bed. For Multi-spouting, the spout diameter increases with increasing bed height, 
however the bed becomes unstable when the bed thickness is increased to 200 mm. For Single¬ 
spouting, the fountain height decreases with increasing bed thickness, but the fountain width and the 
solid volume concentration in the fountain increases. The results of this study provide important 
information on the flow behavior within the multiple-spouted bed and may be helpful for better 
application of this type of spouted bed to the industrial process. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Spouted beds can achieve an intensive gas solid contact and 
can operate stably in a broad range of gas flow rate due to the 
distinguished heat and mass transfer properties and the ability to 
process coarse, sticky, irregularly shaped and heat-sensitive 
materials. In recent years spouted beds have been widely used 
in various industrial and physical operations such as drying of 
biomaterials, foods and pharmaceutical powders (Devahastin 
et al., 1998); coating application on hard and soft capsules 
(Oliveira et al., 2005; Pissinati and Oliveira, 2003); granulation; 
combustion of coals, sawdust and waste biomass (Konduri et al„ 
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1999); coal gasification (Tsuji and Uemaki, 1994); thermal cata¬ 
lytic process (Kechagiopoulos et al., 2009; San Jose et al„ 2009); 
pyrolysis of plastics (Artetxe et al„ 2010; Elordi et al., 2011) and 
flash pyrolysis of biomass (Amutio et al., 2011, 2012). Researchers 
have used and developed various measuring techniques to study 
the flow behavior in the spouted bed, such as non-intrusive 
measuring technique (Mohs et al., 2009), particle tracking 
method, microwave heating and infrared thermal imaging tech¬ 
nology (Zhong et al., 2010) etc. Meanwhile, some experimental 
works have been done to investigate the combustion behavior 
(Pimchuai et al., 2010; Shen et al., 2009) and the gasification 
performance (Spiegl et al., 2010) in the spouted bed reactor. 

In order to satisfy different purposes for using the spouted 
beds and achieve the optimum performance in various industrial 
applications, numerous modifications have been carried out to 
overcome some limitations of the conventional spouted beds. 
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Such as conical spouted bed (Barrozo et al., 2010; Bettega et al., 
2009; San Jose et al., 2008), spout-fluid bed (Littman et al., 2009; 
Patterson et al., 2010), spouted bed with draft tube (Freitas and 
Freire, 2002; Olazar et al., 2012; Shirvanian and Calo, 2004), 
rotating spouted bed (Devahastin et al„ 1999), slot rectangular 
spouted bed (Freitas et al., 2004) and multiple-spouted bed etc. 
However, the scaling-up problem has not been solved well, and 
this limits the application of spouted bed technique in industrial 
processes, especially on large-scale. The multiple-spouted bed 
seems to be a valuable modification in industrial application. 
Several papers (Albina, 2003; Ren et al., 2010; Saidutta and 
Murthy, 2000) about the characteristics and applications of 
multiple-spouted bed have been published. 

To investigate the particle and gas flow characteristics in 
spouted beds, one can choose the experimental approach or 
numerical method. In recent years, numerical simulation with 
computational fluid dynamics becomes very popular to obtain 
detailed information of the flow characteristics in spouted beds 
due to the rapid development of numerical computational tech¬ 
niques and computer performances over past decades. 

The approaches most commonly used in CFD models for the 
simulation of gas solid flow are mainly divided into two groups, 
the Lagrangian approach (discrete models) and the Eulerian 
approach (continuum models). Both approaches consider the fluid 
phase as continuum. The main difference between them is how to 
deal with the dispersed phase. The Lagrangian approach which 
treats particles individually has two branches, the stochastic 
trajectory method and discrete element method. The trajectory 
method is often applied to dilute systems. The discrete element 
method can deal with the dense system with high particle 
concentration and shows a good performance (Link et al., 2009; 
Takeuchi et al., 2004, 2008). However, the computation workload 
increases very quickly with the increase of particle number. 
Consequently, the particle number is too huge to do calculation 
by this model for large scale industrial spouted bed. 

The Eulerian description considers the dispersed phase as a 
continuum. In this approach, the gas phase and solid phase are 
mathematically treated as interpenetrating continuum. Thus, 
supplementary closure equations are required for particle inter¬ 
actions. Subsequently an extension of the classical kinetic gas 
theory called the kinetic theory of granular flow has been 
developed for the closure and is widely accepted (Ding and 
Gidaspow, 1990; Gidaspow, 1994; Lun et al„ 1984). Compared 
to the discrete element method, the Eulerian approach makes the 
scheme of the mathematical solutions relatively easier and 
computational cost less expensive, since the dispersed phase is 
treated with the same discretization and similar numerical 
techniques as those used for the continuous phase, and thus 
greatly benefits the practical applications of this approach. 
Numerical studies (Cunha et al., 2009; da Rosa and Freire, 2009; 
Gryczka et al., 2009; Reuge et al., 2008; Sobieski, 2010) have 
shown the capability of the continuum model incorporated with 
the kinetic theory approach for modeling gas solid flow and have 
demonstrated its superiority for hydrodynamic simulations. 

The two-fluid model based on Eulerian approach assumes that 
both the carrier fluid and the particles comprise two separate, but 
intermixed, continua. For simulating the gas solid flow in spouted 
bed, the motion of particles depends on the proper description of 
the interphase momentum exchange. The interphase momentum 
transfer which has the primary effect on the hydrodynamic 
behavior between the phases is one of the most significant terms 
in the momentum equation, and this momentum exchange is 
represented by drag force. Firstly, the drag force is based on a 
single spherical particle, but people seldom do the simulation of 
just one particle. Secondly, some correlations have been obtained, 
such as Ergun (1952) equation and Wen and Yu (1966) equation. 


However, their validity is limited to the given range of solid 
concentration. In order to overcome this problem and constitute a 
drag equation which can cover a broad range of solid concentra¬ 
tion, Gidaspow combined the Ergun equation and Wen and Yu 
equation, then propose a Gidaspow (1994) model to calculate the 
drag force. This Gidaspow model seems to be more attractive and 
active. It has been widely used and achieved good agreement with 
the experimental findings. 

In this paper, we adopt the two-fluid Eulerian model to 
simulate the complex gas solid two phase flow in a spouted 
bed. The Gidaspow drag model is chosen to describe the inter¬ 
phase momentum exchange. Among the various types of spouted 
beds, we consider a rectangular multiple-spouted bed with a 
central inlet jet and two auxiliary spouting nozzles. The gas solid 
flow patterns in this spouted bed is simulated and compared with 
the recent experimental data. The results show that the overall 
flow behavior and the total bed pressure drop within the spouted 
bed can be predicted well by the aforementioned model. Then, 
the detailed characteristics of the gas solid flow in the multiple- 
spouted bed are examined. The effects of the hydrodynamic 
parameters on the flow patterns in the bed under different 
operational conditions are analyzed, and some interesting and 
useful results are obtained. 


2. Mathematical model 

The two-fluid model based on Eulerian approach is used for 
simulating the gas solid flow in the spouted bed, both phases are 
treated as interpenetrating continua, identified by their volume 
fraction and exchanging properties between the phases. The 
equations of mass conservation and momentum conservation 
are used to describe each of these continua. These governing 
equations are summarily given below. 

2.1. Continuity equations 


The continuity equations for gas and solid phase are 


|(%P g )+Vx ( E gp g u g ) = 0 

(1) 

|(8 s p s ) + Vx( Es p s u s ) = 0 

(2) 


where £ g ,p g ,u g are the volume fraction, density, velocity vector of 
gas phase, respectively, and e s ,p s , u s is the volume fraction, 
density, velocity vector of solid phase, respectively. The additional 
condition for the volume fraction is 

e g +£ s = l (3) 


2.2. Momentum equations 

The gas phase momentum equation has the following form 
| t (w%)+V x (e g PgU g u g ) =-e g Vp+E g p g g+ 

V x Tg + /;(u s -Ug) (4) 

where p stands for the pressure, g stands for the gravitational 
acceleration, x g stands for the gas phase stress tensor, [S denotes 
the drag coefficient. The gas phase is treated as Newtonian fluid 
and the viscous stress tensor is calculated according to Newton’s 
law. 

Tg = -|figPg(V X Ug)I+egPg[Vu g +V T Ug] (5) 


Y. Li et al. / Chemical Engineering Science 80 (2012) 365-379 


367 


where fi g stands for gas phase viscosity, I stands for identity 
matrix. 

The solid phase momentum equation has a similar form to that 
for the gas phase. 

I (w) +v x ^ UsUs ) = -*vp+*p*+ v 

X z s + P(Ug-u s ) (6) 

where t s denotes the solid phase stress tensor, and it will be given 
in the following section. 


2.3. Kinetic theory of granular flow 


For the solid phase, the constitutive equations based on the 
kinetic theory of granular flow are popularly used. This theory is 
derived from the kinetic theory of gas, and a new concept, namely 
granular temperature, is proposed. Compared to the thermody¬ 
namic temperature, the latter represents the fluctuating energy of 
the molecules on the micro-scale; the former expresses the 
kinetic energy of the particle random motion on the macro-scale. 
The meanings of the two concepts are very similar. Thus, the fluid 
dynamic properties of the solid phase, such as the solid pressure 
and viscous force, can be described as a function of the granular 
temperature. The granular temperature 8 is defined by the 
following equation (Gidaspow, 1994). 

e = ^(c 2 ) (7) 

where c is the particle fluctuating velocity, the brackets denote 
ensemble averaging. The transport equation of the granular 
temperature derived from the kinetic theory of granular flow is 

\ [;J( 8 ^)+v x Ow«.0)] = (-p s i+0 

: Vu s + V x (K s V0)-y s +4> s +D gs (8) 

where (-p s I+x s ) : Vu s is the generation of energy by the solid 
stress tensor, p s , is the granular pressure, k s V() is the diffusion of 
energy, k s is the conductivity of the particle fluctuating energy, y s , 
is the collisional dissipation of energy, rp s is the energy exchange 
between the gas and solid phase, and D gs is the rate of energy 
dissipation per unit volume. The conductivity of the particle 
fluctuating energy based on the effective coefficient of restitution 
is calculated by the expression of Gidaspow et al. (1992) as 
follows. 


\50p s d s Vn6 

384g 0 (l+e s ) 


+ 2p s e s 2 d s (l +e s )g 0 


£ 


(9) 


where e s is the coefficient of restitution for particle collisions, d s is 
the particle diameter, and g 0 is the radial distribution function. 
The radial distribution function is used for measuring the prob¬ 
ability of particle collisions, it can be expressed by Ogawa et al. 
(1980) 



( 10 ) 


where s s ,max is the maximum packing volume fraction of particle. 

For granular flow, the particle pressure is calculated by the 
following equation (Lun et al., 1984) which contains the radial 
distribution function. The volume fraction of solid phase is limited 
by its allowed maximum value. When the solid volume fraction 
reaches the packing limit, the particle pressure becomes infinity, 
which implies that the unphysical phenomenon happens. 


( 11 ) 


The solid phase stress tensor can be calculated by the following 
equation 

T 5 = es[(V|ft)(V x u s )jl+e SJ u s [Vu s + V T u s ] (12) 

where A s stands for granular bulk viscosity, p s stands for solid 
phase shear viscosity. The granular bulk viscosity which accounts 
for the resistance to compression and expansion of the granular 
particles is calculated by the following correlation of Lun et al. 
(1984). 

4 C9 \ 1/2 

4= 5 £ sPs d s g 0 (l+e s )(|j (13) 

The shear viscosity of solid phase contains two parts, a 
collisional shear viscosity and a kinetic shear viscosity, as shown 
in the following equation (Gidaspow, 1994). 

Ms = Ms,col +Ms,kin (14) 

The collisional shear viscosity is calculated by the following 
equation (Gidaspow et al., 1992) 

4 /0\ 1/2 

Ms. col = 5 EsPsdsgoO + e s ) (15) 

and the kinetic shear viscosity is expressed by Gidaspow et al. 
(1992) 


Ms, kin 


10p s d s Vnd 
96e s g 0 (l + e s ) 


(16) 


The collisional dissipation of energy represents the rate of 
energy dissipation due to collisions between particles. It is 
represented by the following equation (Lun et al., 1984). 


y s = 


12(l-e 2 )e 2 

dsVn 


PsSo0 3 ' 2 


(17) 


The transfer of the kinetic energy of random fluctuation 
between the gas and solid phase is calculated by the following 
equation (Gidaspow et al., 1992). 

rj> s = —306 (18) 


The last term in the transport equation of the granular tempera¬ 
ture is calculated using the Koch’s expression (Koch, 1990). 


Dgs 


dsP s 

4-Jn8 



(19) 


2.4. The coefficient of interface momentum exchange 

The interface momentum transfer, as mentioned above, is 
represented by the drag force. Since its essential effect on the 
hydrodynamic behavior of gas-solid flow, a lot of work has been 
done on this subject and various drag models have been pro¬ 
posed. As described above, different models including Schiller and 
Naumann drag model, Ergun model, Wen &Yu model and Gidas¬ 
pow model have been put forward. In these models, the Gidaspow 
model has proved to be effective for modeling the interface 
momentum exchange, and has been widely used for simulating 
the solid-fluid systems. Researchers have done some modification 
based on the Gidaspow theory in recent years, one example is the 
Gidaspow model with a switch function (Lathouwers and Bellan, 
2001; Sobieski, 2009; Zhonghua and Mujumdar, 2007). We adopt 
the Gidaspow model in present work and it has the following 
expression (Gidaspow, 1994). 

P= | %(1 d ~ £g) M g | u g-u^|CDo£ g ~ 2 - 65 e g 


p s = B s p s 9+2E s 2 p s () +e s )g Q 0 


>0.8 


(20a) 
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ft =150 ^ % ),gt '4l.75- 
s g d s 


^fig)pg\vg-ns 


where 

Coo = ^(1 +0.15Re 0 687 ) Re < 1000 


Coo = 0.44 Re >1000 


Rc E g p g d s \u g -u s \ 
p e 


e g < 0.8 


(20b) 


(21a) 

(21b) 

( 22 ) 


Table 1 

Parameters and simulation settings in present work. 


Parameters 

Description 

Value 

p s (kg/m 3 ) 

Solid density 

2400 

Pg (kg/m 3 ) 


1.225 

Pg(Pa s) 


1.7894E —5 

d s (mm) 


1.0,1.4,1.8 


Maximum solid volume fraction 

0.59 


Particle restitution coefficient 

0.9 

D 0 (mm) 

Spout nozzle width 

10 

D t (mm) 


300 

H 0 (mm) 

Static bed height 

100,150,200 

H, (mm) 

Bed height 

1200 


3. Numerical procedure 

The experimental results obtained by Ren et al. (2010) are 
adopted to validate the model. The structure of the multiple- 
spout vessel is schematically shown in Fig. 1. The vessel has a 
width of 300 mm. a depth of 30mm and a height of 1200 mm. The 
detailed description of the experimental parameters can be found 
in the paper of Ren et al. (2010). Then, the flow characteristics in 
the multiple-spouted bed under different operation conditions 
have been studied. The parameters used in the present simulation 
are listed in Table 1. We adopt a two-dimensional computational 
grid to simulate the multiple-spouted bed. The effect of the front 
and back walls on the flow behavior in the bed is not considered. 


3.1. Solution method of equations 

The set of governing equations presented in the previous 
sections are solved numerically using a CFD code called MFIX 
(Multiphase Flow with Interphase exchanges) (Syamlal et al.. 
1993) which is developed earlier at the National Energy 

Outlet 


100mm 



Gas inlet 


Fig. 1. Schematic diagram of the multiple-spouted bed. 



x(m) 


Fig. 2. Comparison of the simulated voidages and particle velocities at three 


Technology Laboratory (NETL) of USA and is available on the 
internet at http://www.mfix.org. 

The code uses a finite volume method to solve the governing 
equations. The coupling between the pressure and the velocity is 
obtained using Patankar and Spalding’s SIMPLE algorithm 
(Patankar, 1980). We employed a staggered grid to partition the 
computational domain. The time step is 1x10 4 s. For each 
residual component, a convergence criterion of 10 3 is fixed for 
the relative error between two iterations. All simulations are run 
for 30 s, and data in the last 20 s are used to calculate the time- 
averaged variables. 

The effect of grid partition on the simulation results is 
examined by comparing the simulation results from three differ¬ 
ent grid sizes. As shown in Fig. 2, the predictions from the mid 
grids and fine grids are the same along horizontal direction. Thus, 
we choose the grid partition with the mid grid size for the current 
simulations. 


3.2. Initial and boundary conditions 

Initially, the whole bed is assumed to hold still at atmospheric 
pressure, both of the gas velocity and solid velocity are set to zero. 
The particles are filled in the bed below the given static bed 
height and the volume fraction of particles is limited to the loose 
packing volume fraction. 

For the boundary conditions, the gas is injected into the bed 
from the three spout nozzles, and the gas velocity is specified at 
the spout inlets. The solid velocity is set to zero at the inlets: Zero 
normal gradient condition and the atmospheric pressure are used 
for the outlet at the top of the bed; The no-slip boundary 
condition is used for gas phase at the left and right walls, and 
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the boundary condition of Johnson and Jackson (1987) is used for 
the solid phase. 


4. Results and discussion 

4.1. Validation of the model 

Fig. 3 presents the snapshots of the simulated three typical 
flow patterns, compared with the experimental images photo¬ 
graphed by Ren et al. (2010) using a CCD imaging system. They 
are the internal jet with bubble (IJB), single spouting (SS), and 
multi-spouting (MS), respectively. Reasonable agreements 
between simulations and experiments are found. Fig. 3a shows 
the flow pattern of internal jet with bubble. This is an unstable 
flow pattern, which occurs at high gas spouting flow rates and at 
relatively high static bed heights as in the multiple-spouted coal 
gasifier. In this flow pattern, small bubbles lift off from the jets 
and agglomerate to form a large bubble at the lower part of the 
bed, then the large bubble grows bigger and bigger. When the 
bubble grows big enough to occupy almost the whole width of the 
bed, periodic slugging occurs, which is an awful disaster that 
should be avoided in industrial process. Therefore, the proper 
operation parameters are very important. Simulation results and 
experimental results in Fig. 3b both indicate that the spout is 
stable and non-pulsating, and the movement of particles is 
smooth. The same characteristics as the conventional spouted 
bed can be obviously observed. Fig. 3c also shows the good 
agreement between simulation and experiment. In this case, the 
auxiliary gas spouting velocity is larger than the central gas 
spouting velocity, and two distinct annular fountains can be 
observed, the central fountain is not formed because the annular 
spouts suppress the development of the central spout. 

Fig. 4 compares the simulated and the experimental total bed 
pressure drops versus the central gas spouting velocity at the 
constant auxiliary gas flow rate. Both the simulation and experi¬ 
mental values appear to have the same tendency. The total bed 
pressure drop increases at first, and then decreases gradually after 
it reaches a maximum value. 

From the above comparisons and discussions, it can be seen 
that the simulated results have reasonable agreements with the 
experimental results qualitatively. This indicates that the present 
models are suitable to simulate the gas-solid two phase flow in 
spouted bed and predict a reasonable tendency. 


4.2. Flow behaviors at different central/auxiliary and auxiliary/ 
central spouting gas proportions 

The basic aim of present modeling about the hydrodynamic 
behavior in the multi-spouted bed is to obtain the detailed 
information on flow structure. 


4.2.1. Effect of central spouting gas flow rate 

Fig. 5 shows the voidage distributions within the spouted bed 
for the central injected gas flow rates of O^Qq-O.SQo at the 
constant auxiliary spouting gas flow rate for the particle size of 
1.4 mm, and the constant CL is equal to 1.1 Qm/, where Qm/is the 
minimum fluidizing flow rate. The flow patterns are remarkably 
different for different central spouting gas flow rates. 

The picture in Fig. 5a shows the voidage distributions for the 
central injected gas flow rate of 0.2CL. It can be seen that there is 
not a spout or fountain in the center at this time, because the 
central gas flow rate is very low. The two fountains formed by the 
high auxiliary spouting gases expand through the surface of the 
bed, and throw the particles to the central part of the bed, 
therefore, restricts the growth of the central spouting gas, then 
a dense zone with the particles moving downward is formed and 
leads to a big flow resistance for the gas in the central region. 
Consequently the gas disperses to the bed in all directions as soon 
as it leaves the central spout nozzle, and it is hard to form a jet. In 
this case, particles accelerated by the auxiliary spouting gas travel 
at relatively high speed, while particles in the central region form 
a packed bed and flow downwards at relatively low speed, 
feeding into the spouts along their entire height. 



against central injected gas velocity at the static bed height of H 0 /D m = 1.0. 



Fig. 3. Comparison between experimental and simulated flow patterns under various operational conditions. The left is experimental images, the right is simulated 
snapshots, (a) H„/D 10 = 3.0, Q,/Q m/ = 3.06, Q„/Q m/ = 1.83; (b) H„/D, 0 = 1.5, Q c /Q m/ = 3.6, Q„/Q m/ = 1.65; (c) H 0 /D m = 1.0, Q,/Q m/ = 1.2, Q„/Q m/ = 3.03. 
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Fig. 5. Snapshots of the gas phase voidage distributions with increased central injected gas velocity for 1.4 mm glass beads, (a) Q c /Q a = 0.2\ (b) Q c /Q ( , = 0.3; 
(c) Qc/Qa = 0.4; (d) Qc/Qq = 0.5; (e) CU/Qa = 0.6; (f) Q^/O, = 0.7; (g) QJO* = 0.8. 



As shown in Fig. 5b, when the central spouting gas velocity is a 
little higher than that in Fig. 5a, a void zone is generated only in the 
lower region along the direction of the central spout nozzle. 
However, the central jet is very small and hardly to reach the 
bed surface since the central gas velocity is not high enough, thus 
the turbulence intensity on the upper central region of the bed is 
weak, therefore, the central region of the bed surface is almost flat. 

From Fig. 5c-g, we can see that, at first, the central gas disperses 
to the bed and some pass the packed region, and then enter the 
auxiliary spouts since the gas has the preference to go through the 
small resistance path. As the central spouting gas velocity increases, 
the central spout growth higher and higher. For the effect of the 
gradually strong central spout, the auxiliary spouts tilt slightly 
towards the wall. For the effect of the wall, the fountains are 
squeezed out of shape, and expand to the center. When the central 
spout is high enough to reach the surface of the bed, at last the 


central fountain is formed. However, the central fountain height is 
not very high due to the influence of the auxiliary spouting gas. 

It is the flowing gas in the bed that drives the particle to move, 
so particle performance is different under different operational 
gas conditions. The distribution of particle velocity in the bed 
under various operational conditions is analyzed as follows. 

Profiles of time-averaged particle velocity appear in Figs. 6 and 7 
for the particle sizes of 1.0 mm and 1.8 mm used in this work at 
various bed levels. For each particle diameter, the auxiliary spouting 
gas flow was held constant, and the proportion of central flow was 
varied from 0.20a to 0.80a. The profiles of particle velocity show 
different variation tendencies with the changing gas spouting 
conditions. San Jose et al. (1998) indicated that the particle velocity 
profiles change with the geometry of the spout and, consequently, 
with the operating variables that affect this geometry. Olazar 
et al.(1998) determined that the radial profile of the vertical 
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0 0.05 0.1 0.15 0.2 0.25 0.3 0 0.05 0.1 0.15 0.2 0.25 0.3 

x(m) x(m) 


Fig. 7. Profiles of time-averaged particle velocities with increased central injected gas velocity at various bed levels for 1.8 mm glass beads.(a) Q c /Q 0 = 0.2; 
(b) Qx/Qa = 0.4; (c) Qc/Qa = 0.6 and (d) Qx/Qa = 0.8. 


component of particle velocity changes with the longitudinal 
position. These phenomena also exist in the present spouted bed. 

It is noteworthy that the particles are accelerated more 
intensively by the spouting gas in vertical direction in the jet 
zone and, consequently, a distinct peak value of particle velocity 
is reached close to the spout axis. At the bottom of the bed, the 
velocity of particle is almost zero in the region between the 
spouts and in the region near the wall, because the particles move 
in horizontal direction until they enter the spout regions. 

It can be seen that, when at the lowest ratio of central gas flow 
to auxiliary gas flow as shown in Figs. 6a and 7a, the profiles of 
particle velocity have two maximum values above the auxiliary 
orifices, and then decreases with the increasing distance from the 
spout axis. In the lower part of the bed, the profiles of particle 
velocity are almost flat in the middle region of the bed along 
the horizontal orientation. Near the bed surface, there are two 
regions of negative particle velocity. It is because that, the gas 
conveys a mass of particles to the bed surface due to the high 
velocity of auxiliary gas spouting, and the particles reach their 
maximum value of velocity, then the particles move towards the 
center of the bed. When arriving at the region close to the top of 
the central spout, the particles change their direction and move 
downwards. 

From Figs. 6 and 7 it can be noted that increasing the 
proportion of central gas increases the local particle velocity in 
the central region of the bed due to the higher spouting gas flow 
rate. In the spout region, the particle velocities increase rather 
dramatically with increasing height. In addition, the difference in 
particle velocity between the top and bottom of the bed decreases 
with increasing QJQa. The particles move downwards in the 
region among the spouts and near the wall. In the upper packed 
section of the bed, the particle velocities decrease with decreasing 
height due to cross-flow of solids from the packed region to the 
spouts. In the lower packed section of the bed, the particle 
velocities show the same trend. The cross-flow is greater in the 


lower portion of the bed, probably because of the higher shear at 
the interface between the spouts and the packed region. 

The gas spouting velocity plays a dominant role in the motion 
of particle, and then the movement of particle has an important 
influence on the particle distribution, consequently, while the gas 
spouting velocity is increased, the effect on the particle volume 
concentration distribution becomes pronounced, as presented in 
Figs. 8 and 9. 

At the bottom of the bed, it is observed that the particles in the 
packed region are dense, and the particles in the spout region are 
dilute. In the dense region, the particle volume fraction almost 
approaches 0.59, and the variation of particle volume concentration 
is too small to be observed. In the diluted spout channel, the 
particle volume concentration is very low since the position is 
occupied by the gas. Throughout the whole curve of particle volume 
concentration at lower level, it can be seen that the curve is almost 
flat like a straight line in the packed region, and then decreases very 
sharply to a minimum value close to the axis of the spouts. 

It can be noted that increasing the proportion of central gas, the 
central spout becomes stronger and stronger, and occupies more 
and more space. Finally, it grows almost the same as the other two. 
As shown in Figs. 8 and 9, the particle concentration in the spouts 
increases with increasing height, while the trend is reverse for the 
packed region. The increase of particle concentration in spouts is 
due to entrainment of particles from the packed region and loss of 
gas to the packed region along the height of the bed. 

For all particle sizes, the particle concentration in the central 
region of the bed decreases with increasing Q^IQa- The effect of 
central gas on particle concentration is more pronounced in the 
central region of the bed. The area of packed region decreases 
with the increasing central gas flow rate, since more gas enter 
into the bed, and loosen the particles in the packed region, 
facilitating their entrainment into the spouts. In Figs. 8 and 9, 
one can also observe the gap in particle concentration at the axis 
of the spout between the upper and lower sections of the bed. 

























372 


y. Li et al. / Chemical Engineering Science 80 (2012) 365-379 



Fig. 8. Profiles of time-averaged particle volume concentration with increased central injected gas velocity at various bed levels for 1.0 mm glass beads, (a) Q c /Q 0 = 0.2; 
(b) Qc/Qa = 0.4; (c) Qc/Qa = 0.6 and (d) 0,/Q,, = 0.8. 



This indicates that the momentum of spouting gas is limited, and 
the particles are entrained into the spout mostly in the lower 
section by the high speed gas along the height of the spouts. 


The particle concentration in the packed region varies from 
the bottom to the top of the bed, with the greatest variation 
occurring with the largest particles. Probably because the largest 
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particles is hardest to be entrained by spouting gas in vertical 
direction, meanwhile, the loss of gas to the packed region is 
more than the small particles. By comparing the results of the 
three different particle sizes, it seems that the particle diameters 
dose not have a significant influence on the trend of the 
time-averaged profiles of particle concentration in the range 
investigated. 


4.2.2. Effect of auxiliary spouting gas flow rate 

Fig. 10 shows the voidage distributions within the spouted bed 
for the auxiliary injected gas flow rates of 0.2Qc-0.8Qc at the 
constant central spouting gas flow rate for the particle size of 
1.4 mm, and the constant Qc is equal to 1.1 Qm/. The flow patterns 
are remarkably different for different auxiliary spouting gas 
flow rates. 


As shown in Fig. 10a, when the auxiliary gas flow rate is not 
high enough, the bed resembles a slot-rectangle spouted bed. The 
typical three regions can be observed clearly, the spout, the 
fountain and the packed region. As the auxiliary gas velocity is 
increased, the height of the void zone along the auxiliary nozzles 
increases. The void eventually reaches the bed surface at a certain 
threshold auxiliary gas flow rate and the fountain is formed. The 
threshold gas flow rate is about 0.40c, it is a little lower than that 
in Fig. 5. From Fig. 10, it can be seen that the spouts formed by the 
auxiliary gas tilt towards the axis of the bed, but the fountains are 
pushed aside to the wall under the complex influence of the 
central spouting gas and the bed wall. From Fig. 10, it can also be 
seen that near the top of the spout, there is a dense zone 
surrounding the spout axis, and such dense zone also exists in 
the fountain. This is consistent with the experimental work of 
Grace and Mathur (1978). 



Fig. 10. Snapshots of the gas phase voidage distributions with increased auxiliary injected gas velocity for 1.4 mm glass beads, (a) Q a /Q c = 0.2; (b) Q a /d, = 0.3; 
(c) d,/Qc = 0.4; (d) Qa/Oc = 0.5; (e) O,/0r = 0.6; (f) Qa/Qc = 0.7; (g) 0,/Qr = 0.8. 


a i-5- b 1-5 



Fig. 11. Profiles of time-averaged particle velocities with increased auxiliary injected gas velocity at various bed levels for 1.0 mm glass beads.(a) Q a /Qc = 0.2; 
(b) Qj./Qc = 0.4; (c) Qa/Qc = 0.6 and (d) Qa/Qx = 0.8. 
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According to Fig. 10, we can see that, the height of the central 
fountain increases with the increasing auxiliary spouting gas 
firstly until the auxiliary gas flow rate is more than 0.4Q,-, then 
the height of central fountain decreases with the increasing 
auxiliary gas. It can be interpreted that the auxiliary gas hardly 
move upwards in vertical direction at low gas flow rates, part of 
auxiliary gas is diverted into the central spout for reducing flow 
resistance. The central spout is enhanced and has more power to 
form higher fountain. As the auxiliary gas flow rate is sufficiently 
high to reach the bed surface, the loss of gas to central spout 
gradually decreases, and the fountains formed by the auxiliary gas 
have some influence on the central fountain, then the height of 
the central fountain decreases. 

When the auxiliary gas velocity is low, its influence on the 
central spout is very weak. But as the auxiliary gas velocity 
increases, the influence on the central spout is gradually 
enhanced, and this influence can change the shape of the central 
spout and disturb the flow pattern of the central fountain. 
Compared with Fig. 5, we can see that the influence of auxiliary 
gas on central gas is stronger than the influence of central gas on 
auxiliary gas. 

Profiles of time-averaged particle velocity are presented in 
Figs. 11 and 12 for the particle sizes of 1.0 mm and 1.8 mm used 
in this work at various bed levels. For each particle diameter, the 
central spouting gas flow rate was held constant, and the proportion 
of auxiliary gas flow rate was varied from 0.20c to 0.8Qc. 

The simulation results show that, when the auxiliary gas flow 
rate is very low, the particles move upwards fast in the central 
spout zone, the particle velocity decreases with the horizontal 
distance from the central spout axis, then the particles move 
downwards and reach a negative peak velocity near the interface 
of the spout and the packed region. Finally, the particle velocity is 
almost zero near the wall, they do not move anywhere. Along the 
spout axis, the particle velocity increases firstly, and then 
decreases with the increasing height for the two large particles. 


This characteristic trend is the same as those reported in conical 
spouted beds (He et al., 1994a, b). But the particle velocity 
increases monotonically with height for the smallest particle. 

According to Figs. 11 and 12, it can be seen that the particle 
velocity increases with the increased auxiliary gas velocity in the 
region above the auxiliary nozzles. But, in the central spout, the 
particle velocity increases firstly, then decreases when the aux¬ 
iliary gas flow rate greater than 0.4Q-. As the auxiliary gas flow 
rate increases further, it is clear that the particles move upwards 
fast in the spouts and move downwards in the entire region 
except the spouts. When the auxiliary gas flow rate is increased 
too high to loosen the great mass of the particles in the bed, the 
spouts become unstable and oscillate with time. From Figs. 11 
and 12, we can also see that the particle velocity in the central 
spout is lower than that in the other two spouts when the 
auxiliary gas flow rate is higher than 0.8Q,. It can be interpreted 
as follows: the auxiliary gas accelerates and brings more and 
more particles to the freeboard of the bed by increasing the 
auxiliary gas spouting velocity, then, most of these particles fall 
back to the central region of the bed and move downward. The 
flow resistance in the central region increases, and the velocity of 
up particles is reduced. 

Profiles of time-averaged particle volume concentration are 
shown in Figs. 13 and 14 with increased auxiliary gas flow rate at 
constant central spouting gas flow rate for the particle sizes of 
1.0 mm and 1.8 mm used in this work at various bed levels. 

According to Figs. 13 and 14, the particle concentration 
increases with bed height in the spout region, and decreases from 
the bottom to the top of the bed in the packed region. The present 
simulation results share a similar characteristic with the experi¬ 
mental work on the spout-fluid bed by the previous researchers 
(Pianarosa et al., 2000). 

From Figs. 13 and 14, it can be seen that the particle 
concentration increases with increasing auxiliary spouting gas 
flow rate in the upper part of the central spout. This indicates that 
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Fig. 14. Profiles of time-averaged particle volume concentration with increased auxiliary injected gas velocity at various bed levels for 1.8 mm glass beads, (a) Q„/Q t . = 0.2; 
(b) Q a/Qc = 0.4; (c) Qa/Qc = 0.6 and (d) Q a /Qc = 0.8. 


the auxiliary spouting gas launches a cluster of particles into the 
freeboard and subsequently rain down the central region of the 
bed when the auxiliary gas velocity is sufficiently high. 


It is clearly depicted in the two figures that the profiles of 
particle concentration spread out with increasing bed height. In 
the upper part of the bed, the variation of particle concentration is 
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Fig. 15. Fountain height vs. injected gas flow rate for all three particles, (a) CL, = l.lQ m /, CL = 0-0.9CL; (b) (L = 1.1Q m/ , CL, = 0-0.9Q c . 


less obvious along the bed surface and the fluctuation of the curve 
is moderate, nevertheless they are not flat. 

4.3. Fountain height 

Simulated fountain heights for three different particles are 
presented in Fig. 15. The fountain height is the distance from the 
bottom of the bed to the fountain top. 

Fountain height is a significant characteristic of spouted bed. 
Sobieski (2008) has analyzed sensitivity of Eulerian multiphase 
model for a spouted-bed grain dryer, and the fountain height was 
chosen to be the basic value characterizing the bed. Fie found that, 
besides the inlet gas velocity, the particle diameter was a key 
parameter influencing the fountain height. Even a small modifica¬ 
tion of that value resulted in a significant and clearly noticeable 
change in the fountain height. 

From Fig. 15, it is observed that fountain heights increase with 
increasing spouting gas until the ratio of gas flow rate is greater 
than 0.4, then decrease with increasing spouting gas. This may be 
primarily due to redistribution of fluid between spout and packed 
regions, more fluid flows through the spout when the ratio of gas 
flow rate is lower than 0.4. Because of higher fluid flow through the 
spout, the particle velocity in the spout is also increased, resulting 
in a high fountain height. It is also observed that the fountain 
height fluctuates with the increasing gas flow rate around two 
times the static bed height under the conditions studied. 

At steady spouting, all observed fountain shapes are nearly 
parabolic. But the shape and structure of the fountain has some 
difference from that observed for pure spouting. The fountain in 
pure spouting is extremely dilute, but the fountain in this study is 
observed to have a greater solid concentration. The fountain is 
formed by particle clusters periodically ejected by the bursting 
bubbles. A maximum height is achieved with the expansion of the 
spout when the rising bubbles burst upon the bed surface, and then 
a minimum value is present thereafter. Therefore, the fountain 
height oscillates between a minimum and a maximum value with 
the bubble eruption. Fountain height and shape as well as particle 
velocity and concentration within the fountain has great significance 
to certain applications of spouted bed such as granulators and 
coating units. In view of this, more studies are required to examine 
the effects of operation conditions on fountain characteristics. 

4.4. Spout diameter and shape 

Simulated dimensionless spout diameter as a function of dimen¬ 
sionless height are presented in Fig. 16 for 1.4 mm glass beads under 
varying operating conditions. Dimensionless spout diameter means 
spout diameter normalized by the bed width, and dimensionless 
height means bed height normalized by the static bed height. 
Knowledge of the spout diameter and shape has fundamental 
importance for understanding the hydrodynamics of spouted beds 
and for spouted bed modeling and design. Several researchers (He 
et al., 1998; Mathur and Epstein, 1974; McNab, 1972) have observed 
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Fig. 16. Dimensionless spout diameters as a function of dimensionless height for 
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the spout shape and studied the dependence of spout diameter by 
different experiments with different materials, particle sizes, bed 
dimensions, column geometries and fluid velocities et al. 

From Fig. 16, it can be seen that the simulated diameter of 
spout expands sharply in the region immediately above the inlet 
orifice where the spout diameter is considerably greater than the 
inlet diameter, and then converges slightly near the bed surface. 
The spout shape is related to the interaction between the gas and 
particles at the interface of the spout and packed region. At the 
spout inlet, the gas shears the particle layer at the interface in 
axial orientation and pushes the solids tightly in the horizontal 
direction, then the spout region is enlarged. Ongoing interaction 
between the two phases, the gas carries more and more particles 
as it travels upward in the spout, and then the gas gradually loses 
energy by friction with solids. Therefore, at the upper part of the 
bed, the expanding of the spout begins to decrease, consequently, 
it narrows back. As shown in Fig. 16a, spout diameter decreases 
slightly then remains nearly constant along the height in the 
upper part of the bed. But in Fig. 16b, the decrease of spout 
diameter is considerable. 

The inlet gas spouting velocity has an important effect on the 
spout diameter and shape. The influence of gas flow rate is also 
shown in Fig. 16. It is seen that the spout shape remained the same, 
but the spout diameter increases with increasing gas flow rate until 
the proportion is great than 0.5, then spout diameter decreases 
with increasing gas flow rate. The influence of gas flow rate is more 
obvious in the upper part of the bed than in the lower part. 

4.5. Effect of static bed thickness 

Static bed height is a significant factor affecting the flow 
pattern in the spouted bed, such as spout shape and diameter, 
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Fig. 17. Snapshots of the gas phase voidage distributions with increased static bed thickness for 1.4 mm glass beads, (a) Qc = 0, Q 0 = 1,1 Q m/ ; (b) Qc = l.lQ m /, Q.„ =0. 


fountain height and shape et al. The variation of the voidage in the 
bed with bed height is presented in Fig. 17. 

As shown in Fig. 17a, it appears that the spout diameter 
increases with increasing bed height. However, further increase 
in static bed thickness induces the spouting jet to be submerged by 
the particles below the bed surface, and discharge bubbles peri¬ 
odically. When the static bed thickness is increased to 200 mm, the 
spouting characteristics are quite different from those observed at 
lower bed thickness. Stable spouting becomes unachievable, sub¬ 
sequently the flow pattern in the bed changes to pulsed spouting. 
The fountain collapses, the permanent spout disappears and 
bubbling ensues. The whole bed becomes unsteady and the mixing 
and turbulence of particles in the bed are intensive. 

According to Fig. 17b, with the increasing bed thickness, the 
spout diameter is found to narrow at the middle section but 
expand slightly at the bed surface. As predicted in the equation 
proposed by Grace and Mathur (1978), the fountain height is 
proportional to the square of the particle velocity in the spout at 
the bed surface, therefore higher fountain height will be formed by 
higher gas velocity since higher gas velocity produces higher 
particle velocity. However, the gas velocity in the spout at the 
bed surface decreases with increasing bed height, then for the 
larger bed height, the resulting fountain height is relatively lower 
than for the smaller bed height, but the fountain shape becomes 
wider. In addition, the solid volume concentration in the fountain 
becomes higher and higher with the increasing bed height, 
especially for the central region of the fountain. By Comparing 
Fig. 17a with Fig. 17b, it is found that the effect of static bed 
thickness on multi-spouting is more obvious than single-spouting. 


4.6. Scaling up 

Scale-up of spouted beds is an interesting and important task 
for the implementation of this type of fluid-particle contactors at 
industrial scale. Many scholars have already dealt with this topic. 
Glicksman (1984) proposed a scaling relationship for fluidized bed, 
the controlling non-dimensional parameters were identified as 

w- nr* ft- *£• 5- 

dimensionless particle size distribution, dimensionless bed geometry. 

(23) 

He et al. (1997) modified Glicksman’s (1984) relationship to 
provide a full set of scaling parameters for spout bed scaling up. 
They considered the difference between spouted bed and flui¬ 
dized bed, then added two dimensionless parameters, the internal 
friction angle (<p) and the loose packed voidage (e 0 ) to Glicksman’s 
(1984) relationship by analyzing the force balance for particles 
in the annulus region of spouted bed. The non-dimensional 


parameters were identified as 

W-Mr 1 - Ji’ S' <P’ £ ° 

dimensionless particle size distribution, dimensionless bed geometry. 

(24) 

He et al. (1997) has performed the experimental verifications of 
their scaling parameters, and the modified parameters were valid 
for conventional spouted beds. Although conventional spouted beds 
and multiple-spouted beds share many common features, there are 
also significant differences between them. The flow characteristic in 
the multiple-spouted bed is more complex than conventional 
spouted bed. The multiple spouts interact significantly with each 
other in multiple-spouted bed. Therefore, before the non-dimen¬ 
sional scaling parameters proposed by He et al. (1997) can be used 
to scale up multiple-spouted bed, the detailed validation should be 
conducted. We will do this work in our future study. 


5. Conclusions 

In this paper, the numerical computational technique is used to 
investigate the characteristics of gas solid flow in a multiple-spouted 
bed. Several typical flow patterns are simulated and validated by the 
experimental results carried out by Ren et al. (2010) The simu¬ 
lated results show a reasonable agreement with experimental results. 
Then further investigation has been done when the operational 
conditions are changed. Some significant conclusions are drawn. 

When the central gas flow rate is low, the high auxiliary 
spouting gas completely dominates the flow pattern in the bed 
and restricts the growth of the central spouting gas. With the 
increase of QdQa, the particle velocity increases, but the particle 
volume concentration decreases. The height of central spouting 
jet increases, and the central fountain is formed finally. However, 
the central fountain height is not very high because of the 
influence of the auxiliary spouting gas. The auxiliary gas is 
affected by the gradually strong central spout and deflected 
slightly towards the wall. In addition, the difference in particle 
velocity between the top and bottom of the bed decreases with 
increasing QJQa- On the contrary, the central spouting gas 
controls the whole flow behavior in the bed, and the development 
of the auxiliary gas is suppressed by the central gas. The influence 
of particle size on the flow pattern is unobvious. 

The central fountain height increases with increasing auxiliary 
spouting gas flow rate initially and then decreases. The fountain 
shape is nearly parabolic. The solid concentration in the fountain 
for the current conditions is greater than that for pure spouting. 
The spout diameter expands sharply in the region above the inlet 
orifice and then converges slightly near the bed surface. The spout 
diameter increases with increasing gas flow rate firstly, and then 
decreases. This is similar to the variation of the fountain height. 







378 


V. Li et al. /1 


For the multi-spouting, the spout diameter increases with 
increasing bed height, but the bed becomes unstable when the 
bed thickness is increased to 200 mm, The fountain collapses, the 
spout swings sharply and discharges bubbles. For single-spouting, 
the spout diameter narrows at the middle section but expand 
slightly at the bed surface with the increasing bed thickness. The 
fountain height decreases with the bed thickness, but the solid 
concentration in the fountain increases. 

The hydrodynamic characteristics in multiple-spouted bed with 
mono-dispersed particles has been investigated in this study. How¬ 
ever, the hydrodynamics of spouted beds with poly-dispersed parti¬ 
cles is important in many industrial applications, such as segregation 
and mixing processes. Olazar et al. (1993) and San Jose et al. (1994) 
have done some significant works about the stability and segregation 
of conical spouted beds with binary and tertiary mixtures. For further 
study, the investigation on hydrodynamics in multiple-spouted bed 
when handling solid mixtures will be an interesting topic. 


Nomenclature 

c particle fluctuating velocity, (m/s) 

Cno drag coefficient, dimensionless 

d s particle diameter, (mm) 

D 0 spout nozzle width, (mm) 

D t bed width, (mm) 

Dto bed width of each cell, (mm) 

e s particle-particle restitution coefficient, dimensionless 

g acceleration of gravity, (m/s 2 ) 

g 0 radial distribution function, dimensionless 

Ho static bed height, (mm) 

H t bed height, (mm) 

I identity matrix 

p pressure, (Pa) 

Q c central spouting gas flow rate, (m 3 /s) 

Q a auxiliary spouting gas flow rate, (m 3 /s) 

Q m f minimum fluidizing gas flow rate for a cell, (m 3 /s) 

Re Reynolds number, dimensionless 

u velocity, (m/s) 

u ms spouted-nozzle-based minimum spouting velocity, (m/s) 
U superficial gas velocity, (m/s) 

U a spouted-nozzle-based auxiliary spouting velocity, (m/s) 
U c spouted-nozzle-based central spouting velocity, (m/s) 

Greek letters 

e volume fraction, dimensionless 

e 0 loose packed voidage 

fis.max the maximum of particle volume fraction 

7 S the collisional dissipation term of granular energy 

6 granular temperature, (m 2 /s 2 ) 

2s granular bulk viscosity 

k s granular energy conductivity 

fl the momentum exchange coefficient, (kg/m 3 s) 

p density, (kg/m 3 ) 

t stress tensor 

p viscosity, (Pa s) 

4> s sphericity of particles 

tp internal friction angle of particle phase, (deg) 

Subscripts 

col collision 

g gas phase 


max maximum 
s solid phase 
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